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Abstract 

Magnetic fields of intensities similar to those in our galaxy are also observed in high redshift 
galaxies, where a mean field dynamo would not have had time to produce them. Therefore, a 
primordial origin is indicated. It has been suggested that magnetic fields were created at various 
primordial eras: during inflation, the electroweak phase transition, the quark-hadron phase transi- 
tion (QHPT), during the formation of the first objects, and during reionization. We suggest here 
that the large scale fields ~ fiG, observed in galaxies at both high and low redshifts by Faraday 
Rotation Measurements (FRMs), have their origin in the electromagnetic fluctuations that natu- 
rally occurred in the dense hot plasma that existed just after the QHPT. We evolve the predicted 
fields to the present time. The size of the region containing a coherent magnetic field increased 
due to the fusion of smaller regions. Magnetic fields (MFs) ~ 10/iG over a comoving ~ 1 pc region 
are predicted at redshift z ~ 10. These fields are orders of magnitude greater than those predicted 
in previous scenarios for creating primordial magnetic fields. Line-of-sight average magnetic fields 
(MFs) ~ 10 -2 /iG, valid for FRMs, are obtained over a 1 Mpc comoving region at the redshift z 
~ 10. In the collapse to a galaxy (comoving size ~ 30 kpc) at z ~ 10, the fields are amplified to 
~ 10/iG. This indicates that the MFs created immediately after the QHPT (10~ 4 s), predicted by 
the Fluctuation-Dissipation Theorem, could be the origin of the ~ /j,G fields observed by FRMs 
in galaxies at both high and low redshifts. Our predicted MFs are shown to be consistent with 
present observations. We discuss the possibility that the predicted MFs could cause non-negligible 
deflections of ultra-high energy cosmic rays and help create the observed isotropic distribution of 
their incoming directions. We also discuss the importance of the volume average magnetic field 
predicted by our model in producing the first stars and in reionizing the Universe. 
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I. INTRODUCTION 



The origin of galactic and extragalactic magnetic fields is one of the most challenging 
problems in modern astrophysics [e.g., a a. Magnetic fields on the order of ~ ji G are 
detected in galaxies as well as in clusters of galaxies. It is generally assumed that the coherent 
large scale ~ fiG magnetic fields observed in disk galaxies are amplified and maintained 
by an a — uj dynamo, which continuously generates new fields by the combined action of 
differential rotation (u) and helical turbulence (a). However, the dynamo mechanism needs 
seed magnetic fields and sufficient time in order to amplify them. 

There have been many attempts to explain the origin of seed fields. One of the most pop- 
ular is that they are generated by the Biermann mechanism {3]. It has been suggested that 
;his mechanism acts in diverse astrophysical systems, such as large scale structure formation 
BO , M, cosmological ionizing fronts [7|, and formation of supernova remnants of the first 
stars 8J. Outflows is an additional means of filling protogalaxies with magnetic fields. For 
example, in Section III. G, we discuss outflows of magnetic fields from extragalactic jets, as 
suggested in [9|. 

Another suggestion for the origin of seed fields is that they were created by different 
mechanisms in the very early Universe, before galaxy formation took place. For example, 
such fields may have been created during the quark-hadron phase transition (QHPT), when 
the Universe was at a temperature Tqhpt — W 12 K (Section III. A), during the electroweak 
phase transition (Section III. B), or in the Inflation era (Section III. C). 

One major difficulty with most scenarios for the creation of magnetic fields in the very 
early primordial Universe (-C 1 sec), such as those discussed in §3.1-3.3, is the small coherence 
lengths of the fields at redshifts z < 10. The coherence length is limited by the radius of the 
horizon at the time of the creation of the magnetic field. When expanded to the present time, 
the coherence length is too small to explain the existing observed large coherent magnetic 
fields on the order of the size of galaxies. 

In this paper, we suggest that the observed magnetic fields have their origin in the 
electromagnetic fluctuations in the hot dense plasmas of the very early Universe. This is 
a natural way to create magnetic fields and circumvents the problem of small coherence 
lengths. The Fluctuation-Dissipation-Theorem predicts very large magnetic fields in the 
equilibrium plasma immediately after the QHPT. We evolve these fields to a redshift z ~ 10, 
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when galaxies were beginning to form and find them to be sufficiently strong to explain the 
magnetic field observations in both high and low redshift galaxies. 

We investigate the magnitude of the present magnetic fields in galaxies and the inter- 
galactic medium created by the plasma fluctuations shortly after the quark-hadron-phase 
transition (QHPT), when the plasma properties are well understood. The magnetic fields, 
created by the plasma fluctuations before the QHPT, are poorly understood and we leave 
their evaluation for a future investigation. 

Using the Fluctuation-Dissipation-Theorem (FDT), Opher & Opher [lO, [ll|, Q studied 
the magnetic fluctuations as a function of frequency in the primordial nucleosynthesis era 
and found that they were very large, in particular, at zero frequency. This can be compared 
with the black body prediction which has a zero amplitude magnetic fluctuation at zero 
frequency. 

Tajima et al. [3| suggested that the large magnetic fluctuations predicted by the FDT 
at an early epoch did not dissipate, but continued to exist to the present epoch and now 
contribute to the dominant magnetic field. This scenario is investigated in detail here. Since 
the largest magnetic fluctuations in the plasma occurred shortly after the QHPT, we begin 
our calculations at this epoch. 

Primordial magnetic fields can effect the incoming directions of ultra-high energy cosmic 
rays (UHECRs) above 3 x 10 18 eV. In the last section (Section VI) we discuss the possible 
importance of our predicted primordial magnetic fields on UHECRs. 

We review the observations of astrophysical magnetic fields in Section II Previous sug- 
gestions for creating primordial magnetic fields are given in Section III. The creation of 
magnetic fields in the fluctuations of the hot dense primordial plasma is discussed in Sec- 
tion IV. In Section V , we discuss our model, based on the analysis in Section IV . Our 
conclusions as well as a discussion of our results are presented in Section VI. 



II. OBSERVATIONS OF COSMIC MAGNETIC FIELDS 

The magnetic fields in our Galaxy have been studied by several methods. Measurements 
of the Zeeman effect in the 21 cm radio line in galactic HI regions reveal magnetic fields 
~ 2 — 10/iG. Similar values for the magnetic fields in other galaxies have been obtained from 
Faraday rotation surveys. 
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: Abell clusters have provided information on magnetic 



Observations of a large number o 
ae,d S in clusters of ga^es Q Q M ■ Tne typica! noetic fl e,d strength in the cinster 
is ~ 1 — 10/iG, coherent over 10 — 100 kpc. 

High resolution FRMs of high z quasars allow for the probing of magnetic fields in the 
past. A magnetic field of ~ /i G in a relatively young spiral galaxy at z = 0.395 was measured 
by FRMs from the radio emission of the quasar PKS 229-021, lying behind the galaxy, at z 
= 1.038 Magnetic fields ~ yiiG are also observed in Lya clouds at redshifts z ~ 2.5 [l|. 



III. PREVIOUS SUGGESTIONS FOR CREATING PRIMORDIAL MAGNETIC 
FIELDS 

There have been various scenarios suggested for the source of primordial magnetic fields. 
In this section, we review some of the most important ones. 

A. Magnetic Fields Created at the Quark-Hadron Phase Transition 

In the magnetogenesis scenario at the quark-hadron phase transition (QHPT), proposed 



by Quashnock et al. 18], an electric field was created behind the shock fronts due to the 
expanding bubbles of the phase transition. The baryon asymmetry, which was presumed 
to have already been present, resulted in a positive charge on the baryonic component and 
a negative charge on the leptonic component of the primordial plasma, so that the charge 
neutrality of the Universe was preserved. As a consequence of the difference between the 
equations of state of the baryonic and leptonic fluids, a strong pressure gradient was produced 
by the passage of the shock wave, giving rise to a radial electric field behind the shock front. 
Quashnock et al. Ill estimated the strength of the electric field to be 

eE ~ 15 (-L) (-*-) ( kTQHPT ) f 1 -^-) ™L (I) 
eh ~ 15 U.J VO.l ) \l50MeV )\ I ) cm' {) 

where e is the ratio of the energy density of the two fluids, 5 = (lAp/p), Ap is the pressure 
gradient, and / is the average comoving distance between the nucleation sites. They sug- 
gested that non- negligible fields were produced when shock fronts collided, giving rise to 
turbulence and vorticity on scales of order /. It was found that the magnetic field produced 
on the comoving scale ~ 1 AU has a present magnitude ~ 2 x 10~ 17 G. 
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Cheng & Olinto [19j showed that strong magnetic fields might have been produced 
during the coexistence phase of the QHPT, during which a baryon excess builds up in front 
of the bubble wall as a consequence of the difference of the baryon masses in the quark 
and hadron phases. In this scenario, magnetic fields were generated by the peculiar motion 
of the dipoles, which arose from the convective transfer of the latent heat released by the 
expanding walls. The field created at the QHPT was estimated by Cheng & Olinto [19J to 
be ~ 10 -16 G at the present epoch, on a comoving coherence length ~ lpc. On a comoving 
galactic length scale, they estimated the field to be ~ 10~ 20 G. 

Sigl et al. 20] predicted a present magnetic field ~ 10" 9 G. However, they used very 
special conditions, such as efficient amplification by hydromagnetic instabilities during the 
QHPT. 



B. Magnetic Fields From The Electroweak Phase Transition 

There have been some suggestions made for the origin of primordial magnetic fields 
based on the electroweak phase transition (EWPT). A first order EWPT could possibly 
have generated magnetic fields 20j, |2l| . During the EWPT, the gauge symmetry broke 
down from the electroweak group SU{2) L x U(l)y to the electromagnetic group U{1)em- 
The transition appears to have been weakly first order, or possibly second order, depending 



upon various parameters, such as the mass of the Higgs particle 2l|, |22j. If it were first 
order, the plasma would have supercooled below the electroweak temperature, ~ 100 GeV. 
Bubbles of broken symmetry would have nucleated and expanded, eventually filling the 
Universe. At the time of the EWPT, the typical comoving size of the Hubble radius and the 
temperature were Lfj m 10 cm and T# ~ 100 GeV, respectively. A comoving bubble of size 



Lb — JbLh would have been created with fs — 10~ 3 — 10~ 2 21J. The fluids would have 
become turbulent when two walls collided. Fully developed MHD turbulence would have led 
rapidly to equipartition of the field energy up to the scale of the largest eddies in the fluid, 
assumed to have been comparable to L B . The magnetic field strength at the EWPT would 
have been 

B ~ {Ane) l /\T EW )Tl w {^f-) 2 ~ (7 x 10 21 - 2 x 10 24 )G, (2) 
where e = g^aT^ w /2 ~ 4 x 10 u GeV fm~ 3 is the energy density at the time of the EWPT 



23]. 
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Magnetic fields could also have arisen in cosmological phase transitions even if they were 



of second order 



24l | . In the standard model, the EWPT occurred when the Higgs field 



acquired a vacuum expectation value 77. To estimate the field strength on larger scales, 



Vachaspati [24| assumed that <j) executed a random walk on the vacuum manifold with step 
size £. Over a distance L = N£, where N is a large number, the field changes on the 
average by N 1 ' 2 ^ 1 . On a comoving galactic scale, L = 100 kpc, at the recombination era 
(z~ 1100), Vachaspati 24J found a magnetic field ~ 10 -23 G. 



C. Magnetic Fields Generated During Inflation 

Inflation naturally produced effects on large scales, very much larger than the Hubble 
horizon, due to microphysical processes operating in a causally connected volume before 



inflation (25|. If electromagnetic quantum fluctuations were amplified during inflation, they 
could appear today as large-scale coherent magnetic fields. The main obstacle to the infla- 
tionary scenario is the fact that in a conformally flat metric, such as the Robertson Walker, 
the background gravitational field does not produce relativistic particles if the underlying 



theory is conformally invariant |26|. This is the case for photons, for example, since classi- 
cal electrodynamics is conformally invariant in the limit of vanishing fermion masses (i.e., 
masses much smaller than the inflation energy scale). Several ways of breaking conformal 



invariance have been proposed. Turner & Widrow 25] considered three possibilities: 



1. introducing a gravitational coupling, such as RA^A^ or R^A^A^, where R is the Ricci 
scalar, R^ u the Ricci tensor, and A 11 is the electromagnetic field. These terms break 
gauge invariance and give the photons an effective time-dependent mass. Turner & 
Widrow [25| showed that for some suitable (though theoretically unmotivated) choice 
of parameters, such a mechanism could give rise to galactic magnetic fields, even 
without invoking the galactic dynamo; 

2. introducing terms of the form R^xkF^ F Xk / 'm 2 or RF^F^, where m is some mass 
scale, required by dimensional considerations. Such terms arise due to one loop vacuum 
polarization effects in curved space-time. They can account, however, for only a very 
small primordial magnetic field; and 

3. coupling of the photon to a charged field that is not conformally coupled or anomalous 



coupling to a pseudoscalar field. 

Dolgov & Silk {27] proposed a model invoking a spontaneous breaking of the gauge 
symmetry of electromagnetism, implying non-conservation of the electric charge in the early 
evolution of the Universe. 



D. Generation of the Primordial Magnetic Fields During The Reionization Epoch 

Gnedin, Ferrara & Zweibel [3] investigated the generation of magnetic fields by the 
Biermann battery in cosmological ionization fronts, using simulations of reionization by 
stars in protogalaxies. They considered two mechanisms: 1) the breakout of ionization 
fronts from protogalaxies; and 2) the propagation of ionization fronts through high-density 
neutral filaments. The first mechanism was dominant prior to the overlapping of ionized 
regions (z ~ 7), whereas the second mechanism continued to operate after that epoch as 
well. After overlap, the magnetic field strength at z ~ 5 closely traced the gas density and 
was highly ordered on comoving megaparsec scales. The present mean field strength was 
found to be ~ 10~ 19 G in their simulation. Their results corroborate those of Subramanian 
et al. Q. 

E. Generation of Magnetic Fields Due to Nonminimal Gravitational- 
Electromagnetic Coupling After Recombination 



The generation of magnetic fields by nonminimal coupling was investigated by Opher & 



Wichoski 



29j . From General Relativity, it can be shown that if we have a mass spinning at 



the origin, the space time metric g i is equal to the vector product of the angular momentum 
L and the radial vector r, times 2G/c 3 r 3 , where G is the gravitational constant. Opher &; 
Wichoski 29( suggested that the magnetic field created is proportional to the curl of g^i, 
where the proportionality constant ~ (G) 1 ^ 2 /2c was used, based on the data of the planets 



in our solar system 



3l|. 



Angular momentum in galaxies has 



tidal torques between protogalaxies 



jeen previously suggested to have been created by 



33 



34 



35]. The spin parameter A is defined as 



the ratio of the angular velocity of the protogalaxy to the angular velocity required for the 
protogalaxy to be supported by rotation alone. Numerical simulations find A ~ 0.05, while 
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observations of spiral galaxies show A ~ 0.5. Since A is proportional to the square root of 
the binding energy, it increases by a factor of ten in the formation of a galaxy due to an 
increase of the binding energy by a factor of 100 (i.e., the radius of the protogalaxy decreases 
by a factor of 100). 

In their calculations, Opher & Wichoski 29j investigated models in which the angular 
momentum of a galaxy increased until the decoupling redshift Zd and remained constant 
thereafter. At the decoupling redshift, the spin parameter was A ~ 0.05. They found 
present galactic magnetic fields ~ 0.58 /xG for a decoupling redshift z,i — 100 and noted that 
galactic magnetic fields ~ /iG could be produced by this mechanism without the need for 
dynamo amplification. 



F. Creation of Magnetic Fields From Primordial Supernova Explosions 



Primordial supernova explosions could also be the origin of magnetic fields in the Universe 



0, |36|, |37j] . The scenario investigated was a generic multicycle explosive model, in which a 
Population III object collapsed and then exploded, creating a shock. Matter was swept up by 
the shock, increasing the density by a factor of 4 (for the case of a strong shock). This matter 
was heated to a high temperature, which then cooled down. Eventually spheres of radii of 
approximately half the shell thickness formed and subsequently collapsed into Population III 
stars. They then exploded, starting a new cycle. The supernova shells produced eventually 
coalesced. It was assumed that the gradients of temperature and density in the resultant 
shell were not parallel and that, therefore, a magnetic field was created due to the Biermann 
mechanism. The rate of change of the magnetic field with time is equal to the vector product 
of the density gradient and the temperature gradient times 47r/c£/7ren, where n is the particle 
density and kg is the Boltzmann constant. It was found that this process creates a galactic 
seed magnetic field ~ 10 _16 G, which could be later amplified by a dynamo mechanism. 



G. The Origin of Intergalactic Magnetic Fields Due to Extragalactic Jets 

Jafelice and Opher Q suggested that the large-scale magnetization of the intergalactic 
medium is due to electric current carrying extragalactic jets, generated by active galactic 
nuclei at high z. The action of the Lorentz force on the return current expanded it into the 
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intergalactic medium. Magnetic fields created by these currents were identified as the origin 
of the intergalactic magnetization. They found magnetic fields ~ 10~ 8 G over comoving Mpc 
regions. 

H. Magnetic Field Generation from Cosmological Perturbations 

Another class of magnetic field generation studies are those based on cosmological per- 



turbations. A recent article on this subject is that of Takahashi et al. [38(. They studied the 
evolution of a three component plasma (electron, proton, and photon), taking into account 
cosmological perturbations. The collision term between electrons and photons was evalu- 
ated up to second order and was shown to induce a magnetic field ~ 10" 19 G on a 10 Mpc 
comoving scale at decoupling. 

I. Magnetic Field Generation Due to Primordial Turbulence 

Turbulence has been suggested as the primordial source of magnetic fields. Banerjee & 



Jedamzik 



391 ] has made a detailed study of this scenario. We summarize here their analysis 



and results and compare them with the analysis and results of the present paper. 



It was assumed by Banerjee & Jedamzik [39(] that non-standard out-of-equilibrium 
stochastic magnetic fields were created at high cosmic temperatures T ~ 100 MeV - 100 
GeV, corresponding to quark-hadron or electroweak phase transitions. Their numerical sim- 
ulations were performed using the ZEUS-3D code. Gaussian random fields were used to 
create the non-standard initial turbulent fluctuations. A power-law with distance / was as- 
sumed for the magnetic amplitudes, (B oc l~ n , n = 1-2). The initial stochastic velocity field 
was generated in the same way as the initial magnetic field. A correlation length scale L 
was defined which contains most of the magnetic and fluid kinetic energy. The dissipation 
of the energy into heat occurs via energy cascading from large eddies (~ L) to small eddies 
(~ hiss)- 

Ever since the work of Kolmogorov, it has been known that cascading of energy occurs 
due to eddies on a scale I breaking up into smaller eddies ( ~ 1/2). Typical energy dissipation 
times due to the eddy flows from large to small flows are given by the eddy turnover time on 



the scale L. In the article of Banerjee & Jedamzik 39j, the turnover time on the scale L is 
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comparable to the Hubble time. Thus the turbulent energy introduced in the magnetogenesis 
era is dissipated in one Hubble time. For example, the dissipation time is ~ 10 -4 s for the 
quark-hadron transition. 

The predicted present magnetic field in this turbulent eddy scenario depends on the 
turbulent spectrum assumed at the quark-hadron or electroweak phase transitions. Banerjee 
& Jedamzik 39j found the present magnetic field to be correlated with the comoving 
correlation length L c : B ~ 5 x 10 _15 L C G, where L c is measured in pc. Typically, it was 
found that L c ~ 10 -2 (Eq. (52) in 39]). Thus, the turbulent eddy scenario, with the large 
eddy energy transfers to small scales, where energy dissipation rapidly occurs, typically 
predicts ~ 10~ 16 G on ~ 10~ 2 pc scales. There can be substantial energy transfer to larger 



4Q|. 



scales if the turbulent magnetic field possesses some magnetic helicity 

The above can be compared with the magnetogenesis in the present paper, due to the 
natural fluctuations in thermal equilibrium plasmas. Initially, the magnetic fluctuations 
had an average size A = (77i/3(c/u> p ) [Eq. (15)], where ui p is the plasma frequency and 
c is the velocity of light. They have an average intensity (B 2 )/8n = (T/2)(47r/3)/A~ 3 , 
where T is the temperature. Describing the magnetic fluctuations as dipoles, the magnetic 
field over a distance I due to the randomly oriented magnetic fields follows a power law: 
B = 5(A//) 3//2 . This power law dependence is similar to the power law dependence in the 
turbulent magnetogenesis model, but without a transfer of energy from large to small scales. 
For a thermal equilibrium plasma, the eddy turnover velocity of size I is the thermal rotation 
velocity of the mass of plasma with a diameter /. In the power-law spectrum B = B(\/l) 3 / 2 
with I > X, the eddy turnover time is greater than the Hubble time. There is thus negligible 
energy transfer from the large scale / to the small scale A. On the small scale A, dissipation, 
has already been taken into account by the Fluctuation-Dissipation Theorem. 

The present predicted magnetic field in our magnetogenesis model can be compared with 
the predicted magnetic field of the turbulent eddy model. Whereas, in the turbulent eddy 
model, a present magnetic field ~ 10 -16 G over a comoving correlation length ~ 10~ 2 pc 
is predicted, our model predicts a present magnetic field ~ 10~ 7 G over a comoving length 
~ 1 pc. The predicted magnetic field in our model is , thus, nine orders of magnitude 
greater (over a comoving length two orders of magnitude greater) than that in the turbulent 
magnetogenesis model. This large difference is due to the fast energy transfer from large 
to small dissipation scales in a Hubble time in the turbulent magnetogenesis model, which 
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does not occur in our model. 



IV. CREATION OF MAGNETIC FIELDS DUE TO THE ELECTROMAGNETIC 
FLUCTUATIONS IN HOT DENSE EQUILIBRIUM PRIMORDIAL PLASMAS 



Thermal electromagnetic fluctuations are present in all plasmas, including those in ther- 
mal equilibrium, the level of which is related to the dissipative characteristics of the medium, 
as described by the Fluctuation-Dissipation Theorem (FDT) 4l| [see also Akhiezer et al. 



42 1, Dawson 



43j, Rostoker et al. 44|, Sitenko 45|]. The spectrum of the fluctuations of 



the electric field is given by 



(EiEj)k, 



h 



ku) 



2 e^/ T - 



A 



(3) 



A^,k) = ^f^- 



I +£ij(^,k), 



where £jj(u;,k) is the dielectric tensor of the plasma, u> the frequency, and k is the wave 
number of the fluctuation. From Faraday's law, B = ck/u x E, and setting k = fcx, we find 
for the perpendicular B 2 and B% magnetic fluctuations: 

(B 2 ) kuJ i h c 2 k 2 



8tt 



2 e ^/T-i ^2 



■(A^-A 



33 )i 



and 



(Bl) 



h C 2 k 2 



■(A 



-i 



A: 



2e fa/T-l ,.,2 • 



(5) 



(6) 



87r ■ze" UJ ' J -~ ± u z 

where the subscripts 1, 2, and 3 refer to the x, y, z directions. We then have for the total 
magnetic fluctuations: 



(B 2 ) 



ku 



2u2 



S 7T 



i h c k 
2 e *"/T-i u 2 



(A^ + A^-A 



22 



1* A — 1*> 



(7) 



In order to obtain Aij(u, k) from the equations of motion of the plasma, a multifluid model 
for the plasma is introduced, 

dv, 



dt 



e a E - t] a m a v a , 



(8) 



where a is a particle species label and r] a , the collision frequency of the species. From a 
Fourier transformation of the above equation and rearranging terms, the dielectric tensor 
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can be obtained: 



e y (o;,k) = 5ij - 



CO 



pa 



uj(uj + ir) 



(9) 



where u pa is the plasma frequency of the species a. For, an electron-positron plasma, the 



plasma frequency of the electrons is equal to that of the positrons, u. 



per' 



10 



pe i 



and the 



collision frequencies of the electrons and positrons are equal, rj e 
tensor from Eq. ([9]) then becomes 

,2 



Ve 



(u, k) = Sij — ^—rSij, 

CU{0J + IT}) 



rj. The dielectric 



(10) 



where uJi 



per 



+ u 2 e ^. For electrons, the Coulomb collision frequency is rj e = 2.91 x 



10"VlnAT- 3 / 2 (eV)s- 1 , where n e is the electron density. The collision frequency for the 
case of an electron-proton plasma, which dominates after the primordial nucleosynthesis 
era, is r] p = 4.78 x 10~ 18 n e In AT -3 / 2 (eV)s~ 1 . It describes the binary collisions in a plasma, 
which we assume to be the dominant contribution to i]. We then obtain 



/ 



A; 



1 - 



V 



uj{Lj+iri) 

1 








c 2 k 2 



1 








\ 



c 2 k 2 
,.,2 



(11) 



From Eqs.(7)-(11), the total magnetic field fluctuations as a function of frequency and wave 
number k were found to be 



(B 2 ) k 
8tt 



2hw 



2^2 



k C 



l VUJp (to 2 + 7] 2 )k 4 c 4 + 2u 2 (u 2 p -UJ 2 - T] 2 )k 2 C 2 + \{uj 2 - u 2 ) 2 + V 2 oj 2 ]u' 

(12) 



V. OUR MODEL 



Our model is based on the magnetic fluctuations in the plasma created immediately 
after the QHPT, which are described by the FDT in the previous section. This plasma was 
composed primarily of electrons, photons, neutrinos, muons, baryons and their antiparticles. 
The baryons were essentially stationary and did not contribute to the fluctuations while the 
muons also contributed very little and for a very short time. Since neutrinos are essentially 
massless and act qualitatively like photons, albeit with much smaller cross sections, we 
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assume that they also affect the magnetic fluctuations very little. Therefore we consider 
only an electron-positron-photon plasma before the electron-positron annihilation era and 
an electron-proton plasma thereafter. 

Most of the electromagnetic fluctuations in the primordial plasma that were created 
immediately after the QHPT fall into two broad categories: those with large wavelengths 
(k < u pe /c) at near zero frequency (u <C uj pe ) and those with very small wavelengths (k 

Up e /c) and frequencies greater than u pe . The modes k < u pe /c, denominated "soft" 
or "plastic" photons by Tajima et al. 13j, were significantly modified. It is these plastic 
photons and their magnetic fields in which we are interested. 

From Eq. fll2|) . we obtain the strength of the magnetic field whose wavelengths are larger 
than a size A, 



(5 2 > A /87r = (T/2)(47r/3)A- 3 , (13) 

which decreases rapidly with wavelength. Thus, the magnetic field in Eq.(13) was concen- 
trated near the wavelength A. The spatial size A of the magnetic field fluctuations is related 
to r, the lifetime of the fluctuation, by 

A(r)=2vr-(r /e r) 1 / 2 (14) 

Up 



The average size of the magnetic fluctuations was 

j \[{B 2 ) x /%u}d\ _ 7ir 



A 



Using the model of Tajima et al. 



(B 2 ) x /8-ir}d\ 3 



c/up). (15) 



13j, we assume that a fluctuation predicted by the FDT 
can be described by a bubble of size A. It contains a magnetic dipole whose field intensity 
is given by Eq. ( 1131) . 

The magnetic bubbles were at the temperature of the plasma. We assume that they 
touched each other and coalesced in a time t coa i = \/vb u b, where v^ub was the thermal 
velocity of the bubble. The coalescence time t coa i was found to be much shorter than the 
lifetime r of the bubbles in the primordial Universe, for example, ~ 10 _5 s shortly after t 
~ 10~ 4 s. Before the magnetic fields dissipated, the bubbles coalesced with one another. 
Once a larger bubble was formed, its lifetime, which is proportional to the square of its size, 
was longer. Larger bubbles lived longer and, thus, had more opportunities to collide with 
other bubbles. In this way, a preferential formation of larger bubbles occurred. 

14 



Magnetic field fluctuations are created immediately after the QHPT as predicted by the 
FDT, which we evolve to the recombination era and beyond. Magnetic field fluctuations 
are also predicted to be created by the FDT at the recombination era. Since the created 
magnetic field fluctuations (_B) 2 are proportional to Tn 3//2 , the evolved magnetic fields from 
the QHPT at the recombination era are very much greater than the created magnetic fields 
at the recombination era. The latter source of primordial magnetic fields was thus neglected 



in our investigation. Tajima et al. 



13| previously suggested that the evolved primordial 



fields might continue to exist at the present epoch. No explicit calculation was, however, 
made. Thus, previously it was not known whether these fields would continue to exist or 
not to the present era. We show here that these fields do indeed continue to exit and are not 
destroyed in their evolution by diffusion or reconnection. We also evaluate their structure 
and intensity clS 8b function of redshift. 

We begin our calculations immediately after the QHPT and continue to z ~ 10. Magnetic 
fields were adiabatically amplified at z ~ 10 as the baryon matter collapsed to form galaxies. 

Although the magnetic energy density of neighboring magnetic dipoles is of the same order 
as the energy density of the average magnetic field when they are not at the average distance 
from each other, the magnetic energy density appreciably increases when the neighboring 
dipoles approach each other. Since the field of a dipole is proportional to r~ 3 , where r is the 
distance from the dipole, the magnetic energy density of neighboring dipoles is proportional 
to r~ 6 . Decreasing the separation distance by a factor of 2(4) for example, increases the 
energy density by a factor 64(4096). Thus the magnetic energy density of adjacent magnetic 
bubbles at very short separation distances is very much greater than the average magnetic 
energy density. 

The dipoles tended to align as they interacted due to the intensification of the magnetic 
interaction energy at shorter inter-dipole-distances. The interaction rate of the dipoles 
depended on their thermal velocity. We used as the thermal velocity the velocity of the mass 
of the plasma bubble which is in thermal equilibrium at the temperature of the Universe at 
a given redshift. When the dipoles were oppositely oriented and interacting , two opposing 
processes occurred: alignment and reconnection. As the dipoles approached each other, they 
tended to align in a flip time Tfu p ~ 10~ 5 s shortly after the QHPT at t ~ 10~ 4 s, where Tfu p 
is the time in which a bubble aligns with a neighboring bubble due to the magnetic torque. 

1 /2 

We have Tfu p oc (I/N mag ) , where N mag is the magnetic torque and I is moment of inertia 
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of the bubble. On the other hand, the opposite magnetic fields of the dipoles reconnected in 
a tearing mode time T tear . The shortest T fear is estimated to be f tear = 10°^° t^ 2 t^ 2 , where 



ta = L/va is the Alfven time and tr = 4irL 2 /cr) is the resistive time 46j. The shortest 
tearing time shortly after the QHPT was ~ 10 15 s. Thus Tfu p <C f tear shortly after 10 -4 s 
and remains so for all times of interest. Fig. (1) plots Tfn p , Fig. (2) plots r tear and Fig. (3) 
their ratio, in the time interval ~ 10~ 4 — 10 2 sec. 

The final time plotted in Figs. (l)-(3), ~ 100 s, is the time in which the magnetic field 
in a bubble requires the age of the Universe to diffuse away. Magnetic diffusion, inversely 
proportional to the square of the diameter of the bubble, is only important at early times, 
when the bubbles were small. An initial magnetic field in a bubble diffused away in a time 
Tdif f = 4:7ra L 2 , where L is the diameter of the bubble and a is the electrical conductivity 
47]. 

In the hig h temperature regime (T > 1 MeV) we followed the treatment of Ahonen & 



Enqvist [48[ who solved numerically the Boltzmann equation in the early Universe. For 
T < 100 MeV they found for the conductivity a ~ 0.76T. Since immediately after the 
QHPT the temperature of the Universe was ~ 100 MeV, we used a ~ 0.76T for T > 1 MeV. 
At temperatures T < 1 MeV the conductivity can be approximated as 

a In A \Tim e J 

where A = (1 / 67V 1 / 2 ) (I / a 1 / 2 ) (m 3 Jn e y/ 2 (T/m e ), and a,m e , and n P is the fine structure 



constant, the electron mass, and the electron density, respectively 



49). For L ~ 1 A.U. 



Tdiff is equal to the age of the Universe 47j. In our model the bubbles reached a size ~ 1 
A.U. in a time ~ 100 s. In Figs. (l)-(3) Tfu p and Tf ear are thus plotted from the time of the 
QHPT (~ 10~ 4 s) to ~ 100 s. 

The magnetic field in a bubble would dissipate before coalescence of the bubble occurred 
if the magnetic diffusion time was smaller than the coalescence time. In Fig. (4) we plot 
the ratio of the coalescence time r coa i to the diffusion time Tdiff- It can be seen in Fig. (4) 
that this ratio remains very much less than unity at early times. 

At late times, when the magnetic field flip time (i.e., the time for adjacent dipoles to align) 
was greater than the Hubble time, the magnetic dipoles remained random. The transition 
redshift, when random fields began to exist, was z ~ 10 8 . At this epoch, the comoving size 



1(3 



of the bubbles was ~ 1 pc. In order to explain galactic magnetic fields, we need to evaluate 
the field over the comoving scale of a protogalaxy, ~ IMpc, which eventually collapsed to 
the comoving scale of a galaxy, ~ 30 kpc. 

The magnetic field in a bubble decreased adiabatically as the Universe expanded. Since 
magnetic flux is conserved, we have 

B-h, (IT) 
cr 

where a is the cosmic scale factor. A ACDM model was used to evolve a as function of time, 
with fi m = 0.3, Q A = 0.7, and a Hubble constant h = H/100 km s^Mpc' 1 = 0.72. 

In Figs. (5) and (6), we show the evolution of the size of the bubbles as a function of 
time, from immediately after the QHPT at 10 -4 s to a redshift z ~ 10 at a time ~ 10 16 sec. 
Initially, the size of the bubbles increased rapidly, as shown in Fig. (5). From Fig. (5), we 
observe that the physical size of a bubble increased from 10 cm at t « 10" 4 to 1 cm m a 
time 10~ 7 sec. It continued to increase at this rate until it reached a size ~ 10 7 cm. The 
growth rate then decreased, as shown in Fig. (6). At the redshift z ~ 10 8 (t ~ 3000s), the 
physical size of the bubble was ~ 10 10 cm (i.e., a comoving size ~ 1 pc). 

The manner in which we extrapolated the field amplitude to cosmological scales followed 
the phenomenological analysis of random distributions of size L, in the review article of 
Grasso & Rubinstein [4^] . Their generic average magnetic field over a distance D at a time 
t is proportional to (L/D) p , where p = 1/2, 1 or 3/2, 

(B(L,t)) rms = B (^-) 2 (%X. (18) 



a{t)J \D / 

If we are interested in the volume average magnetic field of a random distribution of dipoles 
in a sphere of diameter D, and each dipole is in a cell of diameter L, the average magnetic 



field is proportional to (L/D) 3 ^ 2 and p = 3/2 in Eq. (fT8l) 231 ] . If, however, we are interested 
in a line-of-sight average magnetic field felt by a cosmic ray particle or a photon (e.g. in 
Faraday Rotation Measurements) the average magnetic field is proportional to (L/D) 1 ^ 2 , 
and p = 1/2 in Eq.fllgjl. 

Non-negligible volume average magnetic field can be important in the formation of the 
first stars and in reionizing the Universe. The formation of the first objects marks the 
transformation of the Universe from its smooth initial state to its clumpy current state. In 
popular cosmological models, the first sources of light began to form at a redshift z ~ 30 
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and reionized most of the hydrogen in the Universe by z ~ 7 50J. In general, it is found 



difficult to reionize the Universe with a standard Salpeter initial mass function 
stellar sources formed by a standard fluctuation dark matter spectrum 



51 



52 



'or the first 



53 



54 



55]. 



Primordial magnetic fields produce additional fluctuations of baryons by the Lorentz force 

scales where the entanglements of 



561 ] . The magnetic tension is more effective on smal 



magnetic fields are larger. Tashiro & Sugiyama 56| found that ionizing photons from 



Population III stars formed in dark halos could easily have reionized the Universe by z 
~ 10 — 20 if the present intensity of the primordial magnetic field is B ~ nG on a comoving 
scale ~ 0.1 Mpc. The relevant Lorentz force causing the collapse of baryon matter is 
proportional to 



(V x B (x)) x B (x) 



D 2 



F. 



(19) 



Thus Tashiro & Sugiyama 56| found that a value F ~ 10~ 28 G 2 /pc 2 is important in 
forming the first objects. In our model a present volume average magnetic field over a 
comoving scale D is B ~ 0.1/zG (D(pc))~ 3 ^ 2 . We thus have F = [10~ 14 / D 2 ]G 2 /pc 2 and 
obtain a Tashiro & Sugiyama 56| value of F with D ~ kpc. We thus find that a D ~ 
kpc comoving region in our model produces a Lorentz force which could be important in 
forming the first stellar sources and in reionizing the Universe. This length is larger than 
the magnetic Jeans length and the cut off length due to direct cascade. Their respective 



wave numbers, given by Tashiro & Sugiyama [56(, are kuj ~ 32Mpc 1 B 1 (nG). and 
k c ~ 102 Mpc- 1 Bq 1 (nG). Putting our volume average magnetic field Bq over D ~ 1 kpc 
into these expressions we obtain kuj ~ 10kpc~ l and k c ~ 34kpc~ l . It is to be noted that a 
sphere of comoving diameter ~ 1 kpc contains a mass ~ 10 3 M Q for a reduced matter density 
Q m ~ 0.3 and Hubble parameter h ~ 0.72. 

A detailed discussion on average procedures of tangled magnetic fields can be found in 
Hindmarsh & Everett 57J. Table (1) shows the growth of the magnetic field in our model 



and the size of the bubbles down to the redshift z ~ 10. The equipartition redshift in Table 
(1) was obtained from the relation (1 + z eq ) « 2.3 x 10 4 f2 m /i 2 [58] . Table (2) shows the 
growth of the line-of-sight average magnetic field over a comoving protogalactic size L ~ 1 
Mpc. 

At z = 10, the intensity of the magnetic field in a bubble whose comoving size is ~ 1 pc 
was ~ Taking the line-of-site average over the comoving scale of 1 Mpc (~ 100 kpc at z 
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~ 10), the rms magnetic field at z = 10 was 9 x 10~ 3 fiG. The magnetic field in the bubbles 
as a function of time is shown in Fig. (7). In Fig. (8), the evolution of the line-of-sight 
average and volume average magnetic field of comoving size ~ 1 Mpc is shown as a function 
of time from i~3x 10 3 sec, when random fields began to exist, to z ~ 10. In the collapse 
of the comoving 1 Mpc region at z = 10 to a galaxy (comoving size ~ 30 kpc), the field is 
amplified to ~ 10/iG. This indicates that the magnetic fields created immediately after the 
QHPT could be the origin of the ~ /xG fields observed in galaxies at high and low redshifts. 



VI. CONCLUSIONS AND DISCUSSION 



We showed that the electromagnetic fluctuations in the primordial plasma immediately 
after QHPT constitute a strong candidate for the origin of primordial magnetic fields in 
galaxies and clusters of galaxies. We calculated the magnetic field fluctuations in the plasma 
after this transition and evaluated their evolution with time. Intense magnetic field fluctu- 
ations on the order of 10 16 G existed at t = 10 -4 sec after the QHPT. These fields formed 
a spatial linkage due to the process of successive coalescence. We showed that magnetic 
bubbles created immediately after the transition could survive to z ~ 10 and could explain 
the observed magnetic fields at high and low redshifts determined by Faraday Rotation Mea- 
surements (FRMs). We found: 1) ~ 10/xG magnetic fields (MFs) over a comoving ~ 1 pc 
region at a redshift z ~ 10; 2) Line-of-sight average MFs, important in Faraday Rotation 
Measurements, ~ 10~ 2 /zG over a 1 Mpc comoving region at z ~ 10 which in the collapse 
to a galaxy comoving size ~ 30 kpc, are amplified to ~ 10/xG; and 3) Volume average MFs 
over a comoving 1 kpc region that could be important in forming he first stellar sources and 



in reioninzing the Universe 



56 



59). 



We found that the magnetic fields in the bubbles, created originally at the QHPT, had 
a value ~ 10/iG at the redshift z ~ 10 and a size 0.1 pc (Table 1). At the present time, 
these bubbles have a comoving length ~ 1 pc and a field ~ 0.1/iG. We can compare these 
results with previous calculations of the creation of magnetic fields at the QHPT. Cheng & 



Olinto 



19], for example, found a much smaller magnetic field, ~ 10 10 yuG, over the same 



also found a much smaller 



comoving size with their mechanism. Quashnock et al. 
resultant magnetic field, ~2x 10~ 11 / uG, over a much smaller comoving size, ~ 10~ 5 pc. 
It is to be noted that the origin of primordial magnetic fields suggested here is qualita- 
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tively different from the other previous suggestions discussed in §3.1-§3.7. These previous 
suggestions require special physical initial conditions. Our model, however, does not. The 
magnetic fields in our model arise from the natural fluctuations in the equilibrium plasma 
that existed in the primordial Universe, described by the Fluctuation-Dissipation-Theorem. 
In 83.8, we discussed the model of Takahashi et al. 



which, like our model, is based 



on natural fluctuations that exist in nature. Our model, however, predicts very much larger 



magnetic amplitudes on a comoving protogalactic scale ~ 1 Mpc. Takahashi et al. [38] found 
a magnetic field B ~ 10~ 25 G on a A = 10 Mpc comoving scale. Since their field is oc k 3 P(k), 
where k is the wavenumber {k = 2tt/X), and the fluctuation power spectrum P(k) oc k n with 



38 



prediction is B ~ 10 G at present for 



n ~ —2 for A < 10 Mpc, the Takahashi et al. 
A ~ 1 Mpc. This can be compared with our prediction for the same comoving scale, which 
is many orders of magnitude greater. 

Our predicted magnetic fields are consistent with present observations. Extragalac- 
tic magnetic fields as strong as ~ 1/iG in sheets and filaments in the large scale galaxy 
distribution, such as in the Local Superclusters, are compatible with existing FRMs 



61 



62 



63 



64 



651 ] . These limits are consistent wit 



i our predicted fields. There is 



66] and a few cases of filaments 



mounting evidence from diffuse radio-synchrotron clusters 

that magnetic fields ~ 0.1 — 1.0/iG exist in the low density outskirts of cosmological 
collapsed objects. These fields may have their origin in the primordial magnetic fields that 
we predict. 

In contrast to the previous models suggested in §3, our model predicts relatively intense 
magnetic fields over small regions in the intergalactic medium. This prediction may help 
to solve the long standing problem of ultra-high energy cosmic rays (UHECRs) (> 3 x 10 18 
eV): their extreme isotropy. The UHECRs are extragalactic since their Larmor radii are 



comparable or greater than the size of the Galaxy [69|, 7(3]. There are only a few nearby 
sources that could be the origin of these cosmic rays. However, the observed arrival directions 
of the UHECRs are highly isotropic jrJQ- 

The importance of intergalactic magnetic fields in creating an isotropic distribution of 
UHECRs has been discussed in the literature. However, different articles arrived at opposite 
conclusions. Whereas Farrar & Piran 173 [argue that the magnetic fields created the observed 

751 ] argue that they are unimportant. Medina- Tanco 



74 



isotropic distribution, Dolag et al. 

& Ensslin To| argue that only weak intergalactic magnetic fields making small angular 
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deflections of the UHECRs may be necessary, since the number of UHECRs sources may be 
much larger than those that are presently observed and that it is possible that fossil cocoons, 
so called radio ghosts, contribute to the isotropization of the UHECR arrival directions 76]. 

Primordial magnetic fields have been previously assumed to exist, without an explanation 
for their origin. Dolag et al. for example, assumed the existence of a homogeneous 

primordial magnetic field ~ 10~ 3 /iG at z ~ 20. They made a magnetohydrodynamic simu- 
lation of cosmic structure formation that reproduces the positions of known galaxy clusters 
in the Local Universe. Protons of energy > 4 x 10 19 eV were found to have deflections, 
which do not exceed a few degrees over most of the sky, up to a propagation distance of 
~ 500 Mpc. It is difficult to explain, however, an isotropic distribution of UHECRs with 
their analysis. 

Relatively intense magnetic fields have been predicted to exist in filaments in the inter- 
galactic medium. Such a filament might exist between us and the powerful radio galaxy, Cen 
A. For example Farrar & Piran [3j suggested that Cen A, at a distance of 3.4 Mpc, could 
be the source of most UHECRs observed. The extragalactic magnetic field was estimated to 
be ~ 0.3/iG. They argue that this scenario can account for the spectrum of UHECRs down 
to m 10 18 ' 7 eV, including its isotropy and spectral smoothness. 

If our predicted magnetic fields are not spread uniformly over space but, as expected, 
are concentrated into the web of filaments, predicted by numerical simulations, appreciable 
deflections of UHECRs propagating along the filaments could occur. The deflection in a 
distance D of a UHECR with energy E = E20 X 10 20 eV by magnetic bubbles of size A and 
magnetic field B is 

6(9) ~ 0.5° [D(Mpc)X(Mpc)} 1/2 B(nG)/E 20 ( 20) 

[73]]. We predicted magnetic bubbles with B ~ 10/iG and A ~ 0.1 pc at z ~ 10. Let us 
assume: that A increased with the cosmic scale factor and that A ~ 1 pc at z ~ 10; and that 
the magnetic fields, trapped in the filaments, decreased slightly to B ~ 1 — 10/uG. From Eq. 
(18), with a distance D ~ 100 Mpc, we obtain 6(9) ~ 5° — 50°/-E , 2o- Appreciable deflections 
could, thus, occur along filaments. 
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FIG. 1: Evolution of the flip time Tfu p (s) of the bubbles as a function of time t(s). 
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FIG. 2: Evolution of the tearing time Tt ear (s) of the bubbles as a function of time t(s). 
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FIG. 3: Ratio of the flip time Tfu p of the bubbles to the tearing time Ttear as a function of time 
t(s). 
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FIG. 4: Ratio of the coalescence time T coa i of the bubbles to the diffusion time Tdiff as a function 
of time t(s). 



25 



2 4 6 8 10 
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FIG. 5: Initial evolution of the physical size of the magnetic bubbles, created immediately after 
the QHPT, as a function of time, t = t Q + At, for t = 10~ 4 s, and < At(10~ 8 s) < 10. 



26 




FIG. 6: Evolution of the physical size of the magnetic bubbles as a function of time from t ~ 0.1s. 
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FIG. 7: Evolution of the magnetic field B{n G) in the bubbles, created immediately after the 
QHPT, as a function of time, t(sec). 
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FIG. 8: Evolution of the line-of-sight average magnetic field B(fi G) of comoving size ~ 1 Mpc as 
a function of time t (sec) from t ~ 3 x 10 3 sec, when random fields began to exist, to t ~ 10 16 sec 
(z ~ 10), when galaxies began to form. 
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TABLE I: Size and Strength of Magnetic Fields in Bubbles 



Epoch 


Magnetic Field (//G) 


Redshift 


Time (sec) 


Size (cm' 


Immediately after the QHPT 


1Q 22 


6 x 10 11 


10" 4 


io- 12 


Electron positron annihilation era 


10 18 


10 10 


1 


10 8 


Nucleosynthesis era 


10 15 


10 8 - 10 9 


1 - 500 


10 10 


Equipartition era 


2 x 10 5 


3600 


10 12 


3 x 10 14 


Recombination era 


2 x 10 2 


1100 


8 x 10 12 


10 15 


Galaxy formation era 


9 


~ 10 


10 16 


10 17 
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TABLE II: Line-of-sight Average Magnetic Fields in Protogalaxies of Comoving Size ~ IMpc 



Epoch 


Magnetic Field (//G) 


Redshift 


Time (sec) 


Size (cm 


Beginning of random fields 


9.5 x 10 11 


10 8 


3 x 10 3 


io- 12 


Equipartition era 


10 4 


3600 


10 12 


10 18 


Recombination era 


300 


1100 


8 x 10 12 


4 x 10 22 


Galaxy formation era 


9 x 10~ 3 


~ 10 


10 16 


10 23 
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